Cesaro's Integral Formula for the Bell Numbers (Corrected) 
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In 1885, Cesaro [1] gave the remarkable formula 

N p = — [ e eC ° sScos(sine)) sin(e cos9 sin(sin6') )smp9 d6 
ire J 

where (A 7 p ) p > 1 = (1,2, 5, 15, 52, 203, . . .) are the modern-day Bell numbers. This formula 
was reproduced verbatim in the Editorial Comment on a 1941 Monthly problem [2] (the 
notation N p for Bell number was still in use then). I have not seen it in recent works and, 
while it's not very profound, I think it deserves to be better known. 

Unfortunately, it contains a typographical error: a factor of p\ is omitted. The correct 
formula, with n in place of p and using B n for Bell number, is 

B n = ^ J e eCOS " cos(sine)) sin(e cose sin(sin6'))sin^ d6 n > 1. 



The integrand is the imaginary part of e e sin n9, and so an equivalent formula is 

2n\ ( r ^ 
B n = Im ( / e e sin n6 d6 

ire 
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The formula (1) is quite simple to prove modulo a few standard facts about set par- 
titions. Recall that the Stirling partition number {^} is the number of partitions of 
[n] = {1,2, ... ,n} into k nonempty blocks and the Bell number B n = ^fc=i {fc} counts 
all partitions of [n]. Thus &!{^} counts ordered partitions of [n] into k blocks (the k\ 
factor serves to order the blocks) or, equivalently, counts surjective functions / from [n] 



1 



onto [k] (the jth block is Since the number of unrestricted functions from [n] to 

[j] is j n , a classic application of the inclusion-exclusion principle yields 

k 
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The trig identity underlying Cesaro's formula is nothing more than the orthogonality 
of sines on [0, %]: 



. 7j if m — n, 
sin m9 sin n6 d6 = < ^ 

if m^n 

for m, n nonnegative integers. Using the Taylor expansion e x = J2 m >o an< ^ DeMoivre's 
formula e ld = cos + % sin 0, it follows that 

Im M e ie ' 9 sin n6> d0j = ^| (3) 



for integer j > 0. Just one more identity is needed: 



/ ^ — *j sinng d 
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for integer & > (of course, {"} = for n > and |£| = for k > n > 0). 
Proof of ( ) The binomial theorem implies the left hand side is 



1 ^(-l) fc -''Qlm Q\ jeie smne d6 

-Y(-i) k - j ( k ) J -- 



0) -i=o 

1 (n \ 7T 
(2) n!\A;j2 
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Summing (1) over > yields Cesaro's formula (1). The Bell numbers have many 
other pretty representations, including Dobinski's infinite sum formula [3, p. 210] 
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